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ABSTRACT 

This paper summarizes the development of an approach 
to optimizing the locations for arrays of sensors and actu- 
ators in active noise control systems. A type of directed 
combinatorial search, called Tabu Search, is used to 
select an optimal configuration from a much larger set of 
candidate locations. The benefit of using an optimized 
set is demonstrated. The importance of limiting actuator 
forces to realistic levels when evaluating the cost function 
is discussed. Results of flight testing an optimized sys- 
tem are presented. Although the technique has been 
applied primarily to Active Structural Acoustic Control 
systems, it can be adapted for use in other active noise 
control implementations. 

INTRODUCTION 

If Active Noise Control (ANC) systems are to be used in 
General Aviation (GA) aircraft, they must be cost effec- 
tive. A primary driver in the cost of ANC systems is the 
number of channels needed to achieve a noise reduction 
target. Additional controller channels increase installed 
cost and subsequent lifecycle costs. The ANC systems 
that have been fielded successfully in several commercial 
turboprop aircraft have dozens of channels of sensors 
(microphones) and actuators (loudspeakers) that are dis- 
tributed throughout the cabin [1, 2], An ANC system 
designed for GA aircraft will be expected to have far 
fewer channels and yet achieve similar noise reduction 
performance. An ANC system with optimized sensor and 
actuator locations can have a reduced channel count for 
a specified noise reduction target and thus be cheaper to 
build and maintain. 


structure to reduce the levels of the structural vibrations 
that cause interior noise. NASA's Langley Research 
Center has investigated the use of ASAC systems to 
control aircraft interior noise [5,6,9]. These results indi- 
cate that actuator and sensor placement strongly influ- 
ence ASAC system performance. However, it has also 
been found that pressurization (at altitude) of the aircraft 
cabin can alter actuator transfer functions and thereby 
reduce the performance of optimized configurations. 

The following sections describe the optimization proce- 
dure and its use in developing an ASAC control system 
for the Raytheon 1 900D. 



Figure 1. Instrumentation layout in Composite 
Cylinder 


ARRAY OPTIMIZATION 

The benefit of optimizing the sensor and actuator arrays 
can be seen in the results of an early test performed on 
NASA's Composite Cylinder [6]. In this test 8 PZT actua- 
tors were installed on the interior shell of the cylinder as 
shown in Figure 1 . Transfer functions were obtained of 


A variant of the ANC system is the Active Structural 
Acoustic Control (ASAC) system [3, 4]. The actuators of 
an ASAC system (e.g., shakers) work directly on the 
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the acoustic response to the individual PZT excitations by 
scanning a microphone boom throughout the interior of 
the cylinder. A total of 462 microphone locations were 
sampled. 

The test focussed on 3 frequencies with specific struc- 
tural and acoustic mode characteristics. These were 210 
Hz (strong acoustic mode, strong structural mode), 230 
Hz (weak acoustic mode, strong structural mode) and 
275 Hz (strong acoustic mode, weak structural mode). 

The goal was to derive and test optimum sensor and 
actuator arrays of dimension 8x4, respectively, at each of 
the test frequencies. For a 4 member PZT actuator array, 
this amounted to finding the best 4 out of 8 possible actu- 
ators, using all 462 microphones. The number of evalua- 
tions required to do an exhaustive search is: 


N a ! 

(N.-N )! x N c ! 


= 70 


(1) 


Where N a =8 is the total number of actuators and N c =4 is 
the number of desired control channels. The exhaustive 
search was done and the best and worst actuator arrays 
were obtained. Finding the optimum 8 microphones out 
of a possible 462 through exhaustive search would 
require a great deal more work: 


( 4 f) = 2.4xl0 : ° 

Clearly an exhaustive search was not feasible. 


(2) 


Random Sensor Locations at 210 Hz 


[□ Best Cose* Worst Case 



Noise Reduction (db) 


Figure 2. Strong acoustic mode, strong structural 
mode 

It is possible to get an estimate of the effect of micro- 
phone optimization. For each of the test frequencies 
(210, 230 and 275 Hz) 1000 random sets of 8 micro- 
phones were selected and evaluated with both the best 


and worst case actuator arrays. The resulting noise 
reduction distributions are shown in Figure 2 -> Figure 4. 


Random Locations at 230Hz 


!□ Best Case* Worst Case 



Noise Reduciton (db) 


Figure 3. Weak acoustic mode, strong structural 
mode 


Random Sensor Locations at 275 Hz 


[□ Best Case * Worst Case 



Noise Reduction (db) 


Figure 4. Strong acoustic mode, weak structural 
mode 

These graphs show how controller performance varies 
with actuator and microphone placement for different 
structural/acoustic properties. In all three cases, the 
microphone placement causes over 5 dB of performance 
variation. When a strong structural mode is combined 
with a strong acoustic mode (Figure 2) actuator place- 
ment is not as important as microphone placement. The 
actuators couple well to the structure wherever you put 
them. With a strong structural mode and weak acoustic 


2 







mode (Figure 3), actuator placement begins to make a 
difference of several dB in controller performance. When 
a weak structural mode is driving a strong acoustic mode 
(Figure 4), actuator placement is vital to good controller 
performance. 

TABU SEARCH 

A form of goal directed combinatorial search call tabu 
search[7], was used to find an optimized microphone 
array. A goal directed combinatorial search uses a cost 
function to selectively refine the target subset using sub- 
stitutions from a larger, candidate set. The danger of 
combinatorial search is that it can get stuck in local min- 
ima that are far removed from the optimum. Tabu search 
provides a means of climbing out of the local minima 
while preventing the search from cycling back. The algo- 
rithm keeps a list of already visited states, the tabu list, 
and will not return to those states. 

The tabu search algorithm can be summarized as fol- 
lows: 

1 ) Select an arbitrary state as Current State. 

2) Loop for Wj Iterations 

a) Put Current State on Tabu List 

b) If Cost Function is minimum, save Current 
State. 

c) Evaluate Cost Function of all States in Neigh- 
borhood. 

d) Move to new Current State not on Tabu List 

which either 

i) Reduces the Cost Function the most, OR, 

ii) Increases the Cost Function the least 

COST FUNCTION REFINEMENTS 

Predicted noise reduction at the error microphones was 
used as a cost function during tabu search. An active 
noise control system can be modeled by: 

e = He + P (3) 

Where P is the primary noise source measured at the 
error microphones, c is the control force input, H is the 
transfer function between the actuator and error micro- 
phones, and, e is the resultant error microphone mea- 
surement. The control force which minimizes (3) can be 
computed by [8]: 


Where, *, denotes the complex conjugate transpose and 
c m is the optimized force vector. The associated cost 
function is a measure of the noise reduction at the error 
microphones and is given by: 

cost = A dB = lOlogg (5) 

The problem with this formulation is that unrealistic actu- 
ator forces can be computed resulting in unrealistically 
high predicted noise reductions and a potentially incor- 
rect optimization. 



Figure 5. Boeing’s Fuselage Acoustic Research 
Facility 

This effect of allowing unreasonably high forces during 
optimization can be seen in results from a noise reduc- 
tion test done in Boeing's Fuselage Acoustic Research 
facility (FARF) [9], see Figure 5. 



c m = -(H*H) ‘h*P 


(4) 
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Figure 6. Single actuator noise reduction, force 
constrained vs. unconstrained 



The FARF is a large anechoic room which contains the 
rear section of a DC-9 aircraft minus engines and tail. 
The interior is complete with seats and trim panels. An 
isolated volume containing 3 rows of seats was formed 
by using 2 acoustically treated barriers. The data were 
originally acquired to support broadband noise control 
experiments. A large external loudspeaker array was 
used as the primary source. A total of 18 error micro- 
phones were located at head height, one for each seat 
(1 5) and 3 in the aisle. Sixty four PZT actuators were 
bonded to the aircraft skin in the frame bays on both 
sides of the aircraft. Only one side of the aircraft section 


<=Front Constrained Case 



was exposed to the primary source. 

Using measured actuator transfer functions, single actua- 
tor noise reduction was computed for each of the actua- 
tors for two cases: actuator force unconstrained and 
actuator force constrained to a reasonable limit, see Fig- 
ure 6. The results are ordered according to constrained 
noise reduction. As can be seen, the unconstrained 
noise reductions are much higher for most actuators, and 
not proportional to their constrained values. The effect 
on optimization can be seen in Figure 7 . 


<=Front Unconstrained Case 



<= Front Speaker Side 


Figure 7. Optimized actuator arrays for constrained and unconstrained cases 


Here 14 actuator locations are optimized for both the con- 
strained and unconstrained cases. The unconstrained 
optimization uses weaker actuators to achieve a greater 
predicted noise reduction at the expense of much greater 
forces. The unconstrained case uses 2 actuators to con- 
trol the area directly opposite the primary source, where 
the constrained case requires at least 8. This illustrates 
how unconstrained forces during optimization can lead to 
incorrect placement of the actuator arrays. 

FORCE LIMITS 

The computed control forces can be limited by the addi- 


tion of a penalty factor, R: 

c m = -(H*H + R) -1 H*P (6) 

The particular penalty factor which corresponds to a 
desired force limit can be found in two ways: (1) an con- 
strained minimization procedure [10], or, (2) conservative 
approximation. The constrained minimization procedure 
most closely resembles the noise control algorithm pro- 
cess where control forces are iteratively increased until 
maximum noise reduction or force limits are reached. 
Although accurate, this procedure is time consuming and 
slows optimization trials. 
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Using a method described by Rossetti [11], the penalty 
matrix, R, is assumed to be a uniform scalar multiple of 
the identity matrix: 

R = U (7) 

Scalar, r, is computed by: 


noise reduction, it has been found to maintain relative 
actuator behavior, selecting the same actuators as the 
constrained minimization procedure. After optimization, 
the final actuator array can be resolved using constrained 
minimization for a more accurate prediction of the noise 
reduction potential. 


Vp*P 7 
a a 2 

max c min 

max 


(8) 


Where a max and a min are the maximum and minimum 
singular values, respectively, of transfer function, H, and 
c max ' s the maximum force limit. The uniform penalty can 
be computed once using all the actuators to form H. This 
value has been found in practice to be very conservative. 


Distribution of Forces - 
CONMIN Method 
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Figure 8. Optimized actuator forces using con- 
strained minimization 

An example solution illustrates the differences between 
the two methods of computing the penalties. An array of 
14 actuators was solved for maximum noise reduction 
over 5 frequencies with the constraint of 5 Vrms maxi- 
mum force. The resultant force distributions are shown in 
Figure 8. As expected, all actuators are operating close 
to maximum force. The same array was solved subject to 
a uniform penalty matrix derived using equations (7) and 
(8). The force distributions are shown in Figure 9. The 
uniform penalty causes an unequal force distribution with 
many actuators way below their maximum potential. 

Although the uniform penalty method does not maximize 


Distribution of Forces - 
Uniform Method 



Figure 9. Optimized actuator forces using uni- 
form approximation 

COHERENCE 

For given coherence, y 2 , between the reference signal 
and the primary source, the maximum noise reduction, 
ignoring force limits, is given by 

A dB = 101og(l- Y 2 ) (9) 

To best predict the noise reduction of an actuator set, the 
effect of coherence must be included in the constrained 
solution derived in the previous section 

Consider the primary sound field to be composed of 
coherent and incoherent parts. Assuming a coherence of 
Y 2 . only the coherent part of primary sound field can be 
controlled. The coherent part of the primary sound field 
is P 00 * 1 and is given by: 

P coh = Py (10) 

The associated coherent error signal is e coh . 

e = Hc + P (11) 

The solution for the optimum, constrained force is now: 
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c m h = -(H*H + R)-’H*P' 


( 12 ) 


coh 


Solving (12) for the optimum constrained force, <4° h , and 
substituting into (1 1 ) will produce the residual coherent 
sound field, e^ h . The total sound field is the sum of the 
coherent and incoherent parts. The incoherent sound 
can be estimated by 

(e inc ) 2 = P*P(1-Y 2 ) ( 13 ) 

The predicted noise reduction now becomes 


A<jb - lOlog 


(e 


coh. 2 


, inc. 2\ 

) + (e ) ' 


P*P 


(14) 


HOW OPTIMIZATION AFFECTS THE CONTROL 
ALGORITHM 


Multi-input/multi-output (MIMO) control systems can be 
plagued by interdependencies between the actuators that 
destabilize the control system. One way of countering 
this effect is to transform the control system into an 
orthogonal coordinate system [12]. This is done by per- 
forming a singular value decomposition on the transfer 
function matrix, H, substituting in (3), and rearranging 
terms. 


USV* = svd(H) 

(15) 

e = USV*c + P 

(16) 

U*e = SV*c + U*P 

(17) 

e = Sv + n 

(18) 


Here, £ and II are the principal coordinates of the error 
microphone array and the primary source respectively. 
The control array, c, is transformed into V. The singular 
values of transfer function matrix, H, are in S. The con- 
troller operates in the coordinate space defined by (18), 
using U* and V to translate the microphone readings into, 
and, the actuator settings out of their respective principal 
coordinates. 



Figure 10. Singular values, S, of 32x12 system 

To illustrate how optimization affects the principal compo- 
nents of a control system, consider an example system 
with 32 error microphones and 12 actuators. Typical sin- 
gular values of the 32x12 system are shown in Figure 10. 
The singular values are scalars and are indicative of the 
authority the virtual control channel, v of (18), has over 
the acoustic power in the associated primary source prin- 
cipal component. For best control, the acoustic power 
should be distributed to take advantage of the increased 
authority in the lower virtual channels. 



Figure 11. Unoptimized principal components of 
the primary source, FL 
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Figure 12. Optimized principal components of 
the primary source, II 

If 12 random locations are chosen for the actuators, the 
principal components of the primary source, H have little 
power in the lowest virtual channels (see Figure 1 1 ). The 
first five virtual channels do hold 75% of the acoustic 
power in the primary source, but, the highest concentra- 
tion of primary power, 20%, is in the 5 th virtual channel. 
When the actuator array is optimized, the primary source 
power becomes more concentrated in the lowest virtual 
channels. In Figure 12, 75% of the acoustic power is 
concentrated in the first 2 virtual channels where there is 
the most control authority as can be seen by comparing 
the distributions in Figure 10 and Figure 12 More noise 
control can be obtained with less power 

RAYTHEON 1900D FLIGHT TEST 

The Raytheon/Beech 1990D, Figure 13, is one of the 
most widely used turboprop airliners in the industry. The 
aircraft can carry 19 passengers 2900 km at a maximum 
cruise speed of 533 kph. The 1900D has a 4 blade pro- 
peller with a blade passage frequency (bpf) of -103 Hz. 
The twin engines are phase locked through a synchro- 
phaser at the shaft speed of 25.8 revolutions per second. 



Figure 13. Raytheon 1900D 


TEST CONFIGURATION 

A block diagram of the system is shown in Figure 13. 
The control system used 32 microphones and 21 actua- 
tor-pairs. The controller, conditioners, amplifiers and dig- 
ital tape recorder were arranged in 2 racks. The trim 
panels and seats were not installed in the aircraft. 

Controller 

The controller consisted of a rack-mounted computer 
with digital signal processor (DSP), I/O and synchro- 
phase interface. The DSP board held two TMS320C40 
processors. The synchrophase interface converted the 
propeller shaft synchrophase signal into a TTL compati- 
ble signal which was then routed into a DSP interrupt. In 
the DSP, an internal timer was synchronized to the inter- 
rupt signal by a software phase-locked loop algorithm. 
The PLL set the internal timer to operate at a multiple of 
48 times the interrupt rate to establish a sampling rate of 
approximately 1238 Hz that would be directly propor- 
tional and locked to the propeller shaft speed (-25.8 rps). 
Filters on-board the I/O boards were set to 723 Hz and 
provided 18 dB roll-off per octave. 
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ANC Block Diagram 


From Synchro Interface 



Figure 14. Block diagram of noise control system 
Control Algorithm 

A principal component least mean squares (PC-LMS) 
algorithm was used as the adaptive control algorithm for 
these flight tests. This algorithm is a transform domain 
version of the multi-channel filtered-x LMS algorithm 
[8,12], and is described in detail elsewhere [13]. In PC- 
LMS the controller parameters (filter weights) are 
adapted in the transformed coordinate system of (18), 
that decouples the feedforward control system at a single 
frequency. Each virtual control channel is independent of 
every other virtual channel. By decoupling the control 
channels, convergence rates and control effort penalties 
can be set for each virtual channel independently. In 
contrast, the filter weights for the filtered-x algorithm are 
adapted in a coordinate system defined by the control 
actuators, which are not usually independent of one 
another and can often show high degrees of inter-chan- 
nel coupling when many actuators are used. 



Figure 15. Error microphone locations 
Acquisition 

The microphones (1/4” electret condenser) were 
clamped to the ring frames and protruded into the cabin 
about 20 cm from the skin. The microphones were uni- 
formly distributed, 4 mics on a ring frame (as shown in 
Figure 15) with the lower and upper microphones roughly 
corresponding to seated and standing head heights, 
respectively. The 8 frames closest to seat locations were 
instrumented. The microphones were connected to ICP 
signal conditioners. 

Actuators 

Motran type IFX 15-100 inertial actuators were installed 
in pairs on the aircraft ring frame. A picture of the IFX 1 5 
is shown in Figure 1 6 and specifications for the actuator 
are summarized in Table 1 . The compact size and high 
force were achieved by the use of Tungsten for the mass. 
The actuator resonant frequency (95 hz) was tuned to be 
just below the 1900D blade pass frequency (103 Hz) to 
avoid the steep phase change that occurs around reso- 
nance. The coil resistance (7.5 £2) was chosen to be 
compatible with the Rane MA 6S multi-channel audio 
amplifiers that were used to power the actuators. 
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Multi-Frequency Performance 



Figure 16. Motran EFX 15-100 



Number of Actuators 


Figure 17. Noise reduction vs. no. of actuators 


Peak Force 

75 N (17 Ibf) 


@ 103 Hz 

Power 

12 W 

Resistance 

7.5 £2 (DC) 

Resonant Freq. 

95 Hz 

Weight 

245 gm (0.5 lb) 

Dimensions 

64x25x36 mm 


(2.5x1 xl .4 in) 


Table 1. Motran IFX 15-100 specifications 


ACTUATOR ARRAY OPTIMIZATION 

Obtaining an accurate characterization of the acoustic 
response due to each actuator is important to the optimi- 
zation process. An actuator must be placed at each can- 
didate location and the acoustic response sampled with 
the microphone array. To avoid the time consuming pro- 
cess of installing actuators at each location, the actuator 
was clamped temporarily in position. The clamped actu- 
ator was able to be fitted to 82 locations on the 1 900D 
frame. Transfer functions were obtained at the bpf and 4 
higher harmonics. All transfer functions were obtained 
during ground tests with the cabin unpressurized. 


Figure 17 shows predicted noise reduction for optimized 
actuator sets of 1 to 24 actuator-pairs. It can be seen 
from the figure that the noise reduction obtained per 
added actuator-pair decreases after 1 2 actuators. Over 
12 dB of attenuation is achieved with 12 actuator pairs 
and just over 14 dB at 24 actuator-pairs. The 12 channel 
system delivers 85% of the noise reduction with 50% of 
the actuators and may be considered a preferred design. 
However, for the purposes of this test all available actua- 
tors (21) were used. 



Figure 18. Optimized locations of 21 actuators 
Figure 18 shows all 21 locations in a view where the 
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frames are unwrapped with bay #1 on the port or left 
side. The 1 st frame is closest to the cockpit, just behind 
the door, and is placed at the top of the figure. Each 
block represents a single possible frame location 
between 2 longitudinal stiffeners. The optimization pre- 
dicted an overall noise reduction of 12.9 dB with reduc- 
tions of 13.5 dB, 8.6 dB and 6.3 dB in the 1 st , 2 nd and 3 rd 
harmonics respectively. 


Harmonic 

Predicted 

Obtained 

1 st 

13.5 dB 

9.5 dB 

2 nd 

8.6 dB 

3.3 dB 

3 rd 

6.3 dB 

1.5 dB 

Overall 

12.9 dB 

8 dB 


Table 2. Noise reduction, predicted vs. obtained 


during the test that a great deal more control effort was 
needed to obtain a corresponding amount of noise con- 
trol. The reason for this can be seen by comparing the 
primary source principal components used during optimi- 
zation with those obtained using the in-flight transfer 
functions. 



Figure 20. Predicted primary source principal 
components 



Figure 20 shows the primary source principal compo- 
nents of the first harmonic for the 21 actuator set using 
the transfer functions acquired during ground tests in an 
unpressurized cabin. The distribution of acoustic power 
is concentrated in the lower virtual channels as expected 
for an optimized array. However, the in-flight principal 
components have shifted, see Figure 21. Over 50% of 
the acoustic power is concentrated in the 3 rd and 4 th prin- 
cipal components. Good noise control is obtainable, but 
at the cost of more control power. 


Figure 19. Reduction of BPF: -measured; -pre- 
dicted; ... predicted using unity coher- 
ence 


RESULTS 

The principal component noise control system worked 
very well providing excellent stability and control man- 
agement. The noise reduction obtained is summarized 
and compared to predicted noise reduction in Table 2. 
The overall noise reduction figure was calculated using a 
linear weighting. 

The actual noise control system did not meet the pre- 
dicted performance levels. This is not due to a lack of 
coherence as can be seen in Figure 19. It was evident 



Figure 21. In-flight primary source principal 
components 
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More extreme shifts in the principal components were 
observed for the 2 nd and 3 rd harmonics. The end result 
is less power available for noise control at all harmonics, 
so lower noise reduction figures are obtained. It is 
believed that this shift in system dynamics is due to cabin 
pressurization. To produce an optimized actuator array 
for an ASAC control system, it is, therefore, necessary to 
use transfer functions acquired under pressurized condi- 
tions. 

CONCLUSION 

Active Structural Acoustic Control (ASAC) been shown to 
be an effective aircraft interior noise control method. The 
ASAC system achieved good control of the blade pas- 
sage frequency. However, the ASAC system did not 
meet performance expectations, especially at higher har- 
monics. This is believed to be caused by changes in the 
dynamics of the structural acoustics of the aircraft 
brought about by cabin pressurization. 

The principal component control algorithm has been 
demonstrated to achieve significant noise reduction. 
The principal component domain transformation has also 
been shown to be a diagnostic tool useful during actuator 
location optimization and post processing analysis of the 
control system. 

The locations for the actuators were optimized using a 
goal-directed combinatorial search. The method coupled 
actuator force constraints and coherence limits to pro- 
duce a realistic prediction of a candidate system's noise 
reduction. The noise reduction predictions were com- 
puted using actuator transfer functions acquired during 
unpressurized ground tests. The use of the unpressur- 
ized transfer functions is shown to bias the optimization, 
producing a sub-optimal actuator array when the cabin is 
pressurized at cruise altitudes, thereby, reducing the per- 
formance of the control system. 

It is finally concluded that if an ASAC system is to reach 
full potential, then the transfer functions used in the con- 
troller and during actuator location optimization, must 
reflect in-flight, pressurized conditions. 
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